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ABSTRACT 

We discuss the characterization of chaotic behaviours in random maps both in 
terms of the Lyapunov exponent and of the spectral properties of the Perron- Frobenius 
operator. In particular, we study a logistic map where the control parameter is ex- 
tracted at random at each time step by considering finite dimensional approximation 
of the Perron-Frobenius operator. 
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1. INTRODUCTION 

The study of systems with many degrees of freedom is one of the central problems 
in the field of dynamical systems [1]. In many cases, there exists a separation of the 
time scales, i.e. a fast evolution superimposed to a slow one, that allows one to capture 
the dynamics in terms of simple models given by one-dimensional random maps [2]. 
For instance, in the fault dynamics, a fast evolution on small scales coexists with a 
slow one on geological times [3] . 

In this paper we shall consider random maps where the fast evolution is taken 
into account by choosing at random at each time step a particular deterministic 
one dimensional map, while the slow evolution is given by the iteration of the map 
extracted time by time. Following the ideas of Spiegel et al. [2] , we start with a system 
made of two variables Xn (representing the slow degrees of freedom) and yn (the fast 
ones), evolving according to the deterministic rule 

{Vn+i =9i{yn) 
(1) 
Xn+l = 92{Xn,yn) 

If the typical evolution time of y„ is much shorter than the corresponding time for Xn, 
such a model can be approximated in terms of the random map 

Xn+l=f^'-\xn), (2) 

where the integer — l,---,k is e.g. an independent random variable, extracted 
at each time step with probability Pi, • • • ,Pn ^ind Z^^-*, • • • , /^'^^ are appropriate one- 
dimensional maps of the interval [0, 1] into itself. This kind of random maps have 
been used to describe many models such as magnetic dynamos [2], transport in fluid 
[4] , random systems [5] and on-off intermittency (abrupt switching among regular 
and chaotic behaviors) [2] . Moreover, random maps exhibit many interesting features 
such as scaling laws for the probability of the laminar phases [2] . 

2 



The characterization of the behaviour generated by random maps requires the 
extension of the tools used for the description of deterministic chaos from a physical 
point of view, and the extension of the rigorous results obtained for expanding maps 
from a mathematical one. In our opinion, both these problems are technically non- 
trivial, and worth to be analyzed because of the physical relevance of the phenomena 
involved. 

The outline of the paper is the following. In sect. II we introduce the basic 
concepts which enter in the characterization of random maps. In particular, we discuss 
the average probability distribution which is the analogue of the invariant measure for 
deterministic dynamical systems, and the so-called snapshot attractor. It is also given 
the definition of the Lyapunov exponent and of complexity for random dynamical 
systems. Sect. Ill is devoted to the Perron-Frobenius operator for random maps. In 
sect. IV we study the case of a random logistic map analyzing the three different 
regimes exhibited by the system for different values of the randomness parameter 
p. We also try to get a deeper understanding of the dynamics by using a finite 
dimensional approximation of the Perron-Frobenius operator. Finally in sect. V we 
draw the conclusions and indicate the possible drawbacks of the Perron-Frobenius 
operator approach to the study of random dynamical systems. 

2. BASIC TOOLS TO STUDY RANDOM MAPS 

Let us briefly summarize the basic concepts needed to characterize random maps. 

1) The average probability distribution. 

After a transient, a deterministic map evolves on the invariant set of the dynamics 
(usually an attractor), where it is possible to define an invariant probability measure. 
We should now consider the probability density obtained as a limit of the histogram 
of points given by the iterations of (2), on the 'average' invariant set /. It is natural 
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to expect that such a distribution might be obtained by means of an average over the 
randomness realizations. Let us recall that for a deterministic map /, the invariant 
probability measure p is the eigenfunction of the Perron-Frobenius operator L related 
to the maximum (in modulus) eigenvalue 71 = 1, i.e. 

p{x)^Lp{x) (3) 

The operator is defined as follows 

L<l>{x)^ fdy6{y-f{x))<P{y)= Yl iTTTVl 

A moment of reflection shows that for a random map it is still possible to flnd the 
average probability density by the straightforward generalization of eq.(3): 

Pav{x) ^ L pav{x) (5) 

where we have introduced the (annealed) average operator 

k 

^=^^3^3 (6) 

and Lj is the Perron-Frobenius operator associated to the deterministic map f^^\ 
2) The snapshot attractor. 

In defining the average probability density, we have considered the long-time 
evolution of a single trajectory under a 'typical' realization of the randomness. One 
can also consider the probability density obtained by starting with a cloud of initial 
conditions x^'^ , Xq^^ , with M very large, that evolve in time under the same 
randomness realization of the dynamics. There exists therefore a instantaneous 
probability density pt{x) obtained by the histogram of the xp'' in the limit M — > oo. 
Such a probability measure is usually indicated as the snapshot attractor [4]. It is 
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generally different of the average probability density, and asymptotically there are the 
following situations: 
(I) pt{x) = 5{x - X*) = pav{x) 
where x* is a stable fixed point; 
(II) pt{x) = 5{x - xl) 

where x^ changes in time and Pavi.^) is not a delta function; 
(III) pt{x) = g{x,t) 

where g{x, t) is a non-trivial function of the space-time. 

The regimes (II) and (III) can be characterized by studying other quantities such as 
Lyapunov exponents, time correlations and complexity measures for random systems, 
see e.g. ref.[6]. 

3) The Lyapunov exponent for random maps 

In noisy systems, the Lyapunov exponent A/ provides the simplest information 
about chaoticity and can be computed considering the separation of two nearby 
trajectories evolving in the same realization of the random process I{t) = ii,i2, it- 
It is possible to introduce the Lyapunov exponent for random maps by considering 
the tangent vector evolution: 



dx 

SO that for almost all initial conditions 



Zn (7) 



Xi = lim ^ln|2;Ar|. (8) 

It is worth to stress that this characterization of the chaoticity of a random dynamical 
system can be misleading [6,7]. A negative value of the Lyapunov exponent computed 
in such a way implies predictability ONLY IF the realization of the randomness is 
known. In other more realistic cases of trajectories initially very close and evolving 



under different realization of the randomness, it can happen that after a certain time 
the two trajectories will be very distant even with a negative Lyapunov exponent A/. 

Both in regimes (/) and (//), the Lyapunov exponent is negative, and in the phase 
(//) it also corresponds to the typical contraction rate of a cloud of points toward the 
'jumping' snapshot attracting point x^, i.e. one has 



On the other hand, we expect A > in the phase (///). 
4) Measures of Complexity 

We pointed out how the characterization of the chaoticity of a random dynamical 
system by the Lyapunov exponent can be misleading. On the other hand it is possible, 
for such a system, to introduce a measure of complexity K which better accounts for 
their chaotic properties [6,7], as 



where hs is the Shannon entropy [8] of the random sequence I{t), A/ is the Lyapunov 
exponent defined above and 9 is the Heaveside step function. The meaning of the 
complexity K is rather clear: ii'log2 is the mean number of bits, for each iteration, 
necessary to specify the sequence xi, ...,Xt with a certain tolerance A. 

We stress again that a negative value of A/ does not implies predictability. To 
illustrate this point, let us calculate K for a system described by a random map which 
exhibits the so-called on- Oj^intermittency [2] (see also sect. IV). In this case one has 
laminar phases, i.e. xt — 0, of average length II and intermittent phases of average 
length //. It is easy to realize that for //<< Zi, 





(10) 



K ^ 




s 



(11) 
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since one has just to compute the contributions of the intermittent bursts whose 
relative weight is ///(// + II) — Ii/Il- 



3. PERRON-FROBENIUS OPERATOR FOR RANDOM MAPS 

A good understanding of random maps can be achieved by studying the spectral 
properties of the Perron-Frobenius operator. In practice, it is convenient to treat the 
problem via finite dimensional approximations. We have that each deterministic map 
f^^\x) {j = 1, - ■ ■ ,k) defines a Perron-Frobenius operator Lj via (4). The evolution 
of the average density is then obtained by the recursive relation 

k 

PU^) = Y.P^^^Pa~v'i^) (12) 

i=i 

Requiring that pav is stationary in time, we see that pav is the eigenfunction of the 
average Perron-Frobenius operator corresponding to the maximum eigenvalue 1, that 
is 

k 

Pav{x) = Lj Pav{x) (13) 

3 = 1 

The Lyapunov exponents can then be obtained by a space average over the attractor 



as 



Xi 

.1 T 
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pav{x)dx (14) 



dx 

Moreover, the spectrum of the Perron-Frobenius operator is directly related to the 
decay rates of the time correlations also called 'resonances' of the dynamical system 
[9]. In particular if we consider a generic smooth observable the time correlation 
at large times decays as 

CA{t) =< A{xt)A{xQ) > - < A >2 - e"^* (15) 

where 7 > 0, exp(— 7) is the modulus of the second eigenvalue of the Perron-Frobenius 
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operator and the average is a time-average, 




(16) 



For random maps 7 is given by the logarithm of the modulus of the second eigenvalue 
of the average Perron-Frobenius operator L. 

Beyond the average operator, it is interesting to consider the product of random 
transfer operators operator 



that is itself a random operator acting in an appropriate functional space. Roughly 
speaking, it can be regarded as the product of infinite dimensional matrices whose 
elements are transition probabilities. It is therefore natural to expect that in the 
limit t ^ 00, {VtP^)'^ for almost all the randomness realizations has 
a non-random spectrum, as it happens for the product of random matrices, where 
the Oseledec theorem [10] holds. This result has been proved in the mathematical 
literature in the case of infinite dimensional operators acting on Banach spaces [11]. 
The operator Vt is interesting since it controls the evolution of the instantaneous 
probability density Pt{x): 



This problem is analogous to that one of the evolution of a tangent vector in dynamical 
systems, see ref.[12]. In that case the jacobian matrix plays the role of the Perron- 
Frobenius operator and the tangent vector the role of the probability density. Note, 
however, that in the usual deterministic dynamical system problems the jacobian 
matrix is not extracted at random, although given by a chaotic dynamics. 

The first eigenvalue of V is equal to unity since pt{x) is a normalized density. The 
second eigenvalue 72 < 1 measures the exponential rate of collapse of a generic density 



Vtih,- ■ ■ ,it) = Lit-i ■■■L 



(17) 



Pt{x) = VtPo{x) 



(18) 
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Pt toward the time-dependent eigenfunction p^. In other terms in an appropriate norm 
one has 

llpt-ptll-e*'"!^^! (19) 

It is worth stressing that in the phase {B) where the attractor is the jumping point 
x^, In 1 72 1 = A measures the contraction rate of a cloud centered around x^. On the 
contrary, in the fully chaotic phase (C), In IA2I has no particular intuitive meaning. 

4. TRANSITION TO CHAOS IN THE RANDOM LOGISTIC MAP 

In this section we analyze a specific example of random map, namely 

xt+i = x{l-x) (20) 

with 

4 with probability p 

(21) 

1/2 with probability 1—p 
The evolution is thus given by the random composition of a contracting and of an 

expanding logistic map. For sake of simplicity the expanding map has the control 

parameter set at the Ulam point (r = 4). Indeed, for p = 1, it is possible to find by a 

well known duality argument 

the invariant probability density, 

p(^)-:^^7=^ (22) 

TT ^x{l - X) 

At varying p, one observes three different regimes: 
(1) p < |. The random map has a stable fixed point in the origin. The average 
probability density Pav{x) — d{x) and the Lyapunov exponent is negative and it 
is trivially given by the expression 

aW-P + A^2).(i_^) (23) 
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where A^^) = 4 and A^^) = I are the Lyapunov exponents of the two logistic maps 
(21). 

(2) ^ < p < Pc = 0.47.... There is a snapshot attractive fixed point and Pav{x) ^ 
on all the interval [0, 1] while pt{x) — d{x — x^). The Lyapunov exponent is still 
negative and at least near the transition the complexity is given by the expression 
already cited in section 2. 

K ~ ^-^hs (24) 

I'L 

where II and Ij are respectively the average lengths of the laminar and intermittent 
phases and the (24) is written in the case // << II (see fig. 1). 

(3) p > Pc. 

The random map is chaotic: the Lyapunov exponent is positive, both Pav{x) 
and pt{x) 7^ on all the interval [0, 1]. 

We have tried to get a deeper understanding on the dynamics by considering finite 
dimensional approximations of the Perron-Frobenius operator. The idea is to treat 
the two Perron- Probenius operators Li and L2 related to the contracting and to the 
expanding map as N x N matrices. In fact, there are two possible strategies. 

The first one is to consider a uniform partition of the interval [0, 1] made of N 
sub- intervals. Then we approximate L by the matrix of the transition probabilities 
between the partition elements. In the limit N ^ 00 such a matrix converges to L 
according the Ulam conjecture [13], that can be proved in some cases [14-15]. 

The second one is to look for the matrix representation of the operator in an 
appropriate base of the functional space of polynomials up to order {N — 1), i.e. the 
so-called Galerkin approximation [16]. 

We have used both the strategies, even if it turns out that the Ulam approximation 
is in practice more convenient. The Galerkin approximation using Legendre 
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polynomials gives unavoidable numerical troubles since the coefficients of the powers 
grow very fast with the order n of the polynomial. 

We performed the numerical treatment of the Perron-Frobenius operator as follows. 
For a given partition of the interval (0, 1) in segments the Perron- Probenius operator 
Lq, {a = 1,2) of the map a is approximated by the N x N matrix Lo,{N) whose 
elements [Lct{N)]ij are given by the probability to perform a transition from the 
interval j to the interval i and L is the matrix 

L{N)^pLi{N) + (1 -p)L2{N). (25) 

Prom (25) one can easily compute Pavi^), and 7. 

Pig. 2 shows A/ vs p. One can observe how the Lyapunov exponent computed 
using the approximated Perron-Probenius operator (25) is in good agreement with 
the exact (numerical) results only in the region (///) and for very small values of p. 
On the contrary the finite dimension approximation of the Perron-Probenius operator 
seems to have (if any) a very slow convergence in the region (//). 

Pig. 3 shows C(t) vs t at different values of p. One roughly has ior p = pc + e 

C(t) ~ e-T'(^)^T-"^^) (26) 

where 7(e) — > for e — > 0. This behaviour recalls a critical phenomenon and it is 
essentially understood [2]. The approximation (25) of the Perron-Probenius operator 
gives very bad results for 7. 

Pig. 4 shows Pav{x) vs X for two different values of p. In the region (II) for 
p = Pc + e one has a power law behaviour 

Pav{x)^X-<'l (27) 

The approximation of the Perron-Frobenius operator is able to give a good 
approximation in the regions (//) and (HI)- 
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5. CONCLUSIONS 



In this paper we have discussed the dynamical characterization of systems whose 
time evolution is described by random maps, with particular emphasis on the notion of 
Lyapunov exponent and of complexity. We have studied the application of the Perron- 
Frobenius operator to the analysis of the dynamical behaviour of random dynamical 
systems. The Perron-Frobenius operator can be analyzed through finite- dimensional 
approximations in order to derive the invariant measure of the system, the Lyapunov 
exponent and the temporal correlation functions. In particular we have followed two 
ways. 

The first way consists in approximating of the Perron-Frobenius operator as an 
N X N matrix whose elements are the transition probabilities between the elements 
of the partition of the interval [0, 1] in N sub-intervals. According to the Ulam 
conjecture such a matrix should converge to the Perron-Frobenius operator in the 
limit A?" — > oo. The second approach is the so-called Galerkin approximation and 
consists in the matrix representation of the operator in an appropriate basis of the 
functional space [15]. We found out that the first strategy gives the better results. 
In particular it provides a good approximation of the average invariant measure pav 
(see fig.s 4a and 4b). As for the value of the Lyapunov exponent A/, the method is 
suitable to describe its behaviour near the transition to chaos (from the region (//) 
to the region (///)) but the results are less good near the first transition between the 
regions (!) and (//). In this case, in fact, the contributions to 



could be large exactly where the Perron-Frobenius method is not able to give an 
accurate approximation for pav Fig. 2 shows that at increasing the order of 
approximations of the Perron-Frobenius operator from a 100 x 100 to a 200 x 200 
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matrix, there is a tendance to obtain better and better approximation of the exact 
value of A/. We thus expect that in the hmit N ^ oo this approach tends to recover 
the real value of A/. As for the temporal correlation functions the method based on 
the approximation of the Perron-Frobenius operator does not work due, probably, to 
a too slow convergency with N. 
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Figure captions 



Fig. 1 xt ys t for the random logistic map described by eq.s (21) with p = 0.35. 

Fig. 2 Lyapunov exponent A/ vs p for the logistic random map described by eq.s(21); 
(a) numerical result. Curves (b) and (c) indicate the theoretical estimation of the 
Lyapunov exponent with the Ulam method with matrices NxN with (b) : N = 200 
and {c):N^ 100. 

Fig. 3 Temporal autocorrelation functions for the random logistic map (21) for three 
different values of p. 

Fig. 4 Histograms of the distribution of positions for the random logistic map. The dot- 
dashed line corresponds to the theoretical estimation and the solid line to the 
numerical result. (4a) p — 0.35, (4b) p — 0.48. 
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